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Abstract. An explicit SUSY SU(12) unification model with three light chiral families is presented which avoids any external 
flavor symmetries. The hierarchy of quark and lepton masses and mixings is explained by higher dimensional Yukawa 
interactions involving Higgs bosons containing SU(5) singlet fields with VEVs appearing at or below the SUSY GUT scale 
of 2 x 10 16 GeV, approximately 50 times smaller than the SU(12) unification scale. The model has been found to be in good 
agreement with the observed quark and lepton masses and mixings, with nearly all prefactors of G(\) in the four Dirac and 
one Majorana fermion mass matrices. 
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MOTIVATION 

Grand Unified Theories (GUTs) of quarks and leptons have played an important role for almost 40 years in theoretical 
attempts to make sense of the apparent group structure and mass spectra of the quark and lepton matter fields and 
coupling strengths of the gauge and Higgs fields observed in Nature. Although many varieties of models have been 
constructed, the most popular unified ones are based on the unitary, orthogonal, or exceptional groups SU(5), SO(10), 
or E6, respectively. But in each of these cases, the chiral irreducible representations (irreps) can uniquely describe only 
one family of quarks and leptons: towit, 10, 5, and 1 for SU(5), 16 for SO(10), and 27 for E6. In order to accommodate 
the three families observed todate, it has been conventional to introduce in addition to one of the above Gf am ;iy groups, 
a Gfi avor symmetry group which also distinguishes the families. While continuous flavor symmetries such as U(l), 
SU(2), SU(3) and their products have been considered in the past, more recently discrete symmetry groups such as 
A4, T' and S4, etc. have been fashionable in the past 10 years [1, 2]. In either case, the GUT model then involves the 
direct product group G family x Gfi avor . 

True family and flavor unification requires the introduction of a higher rank simple group. Some earlier studies 
along this line have been based on SO(18) [3, 4], SU(ll) [5, 6], and SU(9) [7, 8]. More recently, models based on 
SU(7) [9], SU(8) [10], and SU(9) again [11, 12], the latter reference by two of us (RF and TWK), but none have been 
totally satisfactory due to a huge number of unwanted states and/or unsatisfactory mass matrices. Here we describe an 
SU(12) unification model [13] with interesting features that was constructed with the help of a Mathematica computer 
package called LieART written by two of us (RPF and TWK) [14]. This program allows one to compute tensor 
products, branching rules, etc., and perform detailed searches for satisfactory models in a timely fashion. While other 
smaller and larger rank unitary groups were examined, a model based on SU(12) appeared to be the most satisfactory 
minimal one for our purpose. We sketch here the model construction and point out that further details can be found in 
[13]. 

SU(12) UNIFICATION MODEL AND PARTICLE ASSIGNMENTS 

While the three popular GUT groups cited earlier each have just one usefu l chi r al irr e duc ible representation (irrep), 
SU(12) has 11 totally antisymmetric irreps: 12, 66, 220, 495, 792, 924, 792, 495, 220, 66, and 12, of which 10 
are complex (while 924 is real), which allow three SU(5) families to be assigned to different chiral irreps. For this 
purpose, one chooses an anomaly-free set of SU(12) irreps which contains three chiral SU(5) families and pairs of 



fermions which will become massive at the SU(5) scale. One such suitable set consists of 

6(495) +4(792) +4(220) + (66) + 4(12) -> 3(10 + 5 + 1) +238(5 + 5) +211(10 + 10) +484(1) (1) 

where the decomposition to anomaly-free SU(5) states has been indicated. The latter follows from the SU(12) -> 
SU(5) branching rules: 



495 -> 35(5) +21(10) +7(10) +5 + 35(1), 

792 -> 7(5) +21(10) + 35(10) + 35(5) +22(1), 

220 -> 10 + 7(10) +21(5) + 35(1), 

66 -> 10 + 7(5) +21(1), 

12 -> 5 + 7(1) 



(2) 
(3) 
(4) 
(5) 
(6) 



A search through the possible assignments of the three light chiral families to the SU(12) irreps appearing in the 
anomaly-free set of Eq. (1) reveals the following selection for a satisfactory low scale phenomenology: 
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3rd Family : 
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Here the subscripts on the SU(12) irreps refer to the family in question, while the numbers in parentheses are just the 
SU(5) irreps chosen. Note that each SU(5) family multiplet can be uniquely assigned to a different SU(12) multiplet 
in the anomaly-free set according to (1). On the other hand, the remaining SU(5) multiplets are unassigned but form 
conjugate pairs which become massive and decouple at the SU(5) scale and are of no further interest. 



EFFECTIVE THEORY APPROACH AND LEADING ORDER TREE DIAGRAMS 

We start with the SU(12) model sketched above and take it to be supersymmetric. With a 143h adjoint Higgs field 
present, the breaking of SU(12) to SU(5) can occur via SU(12) -» SU(5) x SU(7) x U(l), and in steps down to 
SU(5) via a set of antisymmetric chiral superfield irreps appropriately chosen to preserve supersymmety [15, 16]. 
Unbroken supersymmetry at the SU(5) GUT scale allows us to deal only with tree diagrams in order to generate 
higher dimensional operators, for loop corrections are much suppressed. 

For this purpose, we introduce massive 220 x 220 and 792 x 792 fermion pairs at the SU(12) scale. In addition, 
we introduce (1)66h, (1)66„, and (1)220h, (1)220h conjugate Higgs pairs which acquire SU(5) singlet VEVs at 
the SUSY SU(5) GUT scale. Finally, doublets in (5)924 H and (5)924 H Higgs fields effect the electroweak symmetry 
breaking at the electroweak scale. The list comprises then the following: 

Higgs Bosons Massive Fermions 
(5)924 H (5)924 H 220x220, 

(1)66h (1)66 h 792x792 (8) 

(1)220 H (1)220 H 

(24)143 h 

For each element of the quark and lepton mass matrices, tree diagrams can then be constructed from three-point 
vertices which respect the SU(12) and SU(5) multiplication rules. 

For illustration we present the lowest order tree diagram contributions to the 33 elements for the up and down quark 
mass matrices, taking into account the family assignments in (7). These are listed as U33 and D33, respectively in (9). 



Table 1. Leading order up- and down-type quark diagrams for each mass matrix element. 



Up-Type Quark Mass-Term Diagrams 

( 1 0)220^.(5)92% .(10)22(b 
(10)792 2 .(1)66 h .(10)220x(10)220.(5)924h 
( 1 0)220 3 . (5)924 H . ( 10) 220 x ( 1 0) 220. ( 1 )66 H 
(10)495 1 .(l)220 H .(T0)792x(10)792.(l)66 H 
( 1 0)2203 . (5)924 H .(10) 220 x ( 1 0)220.(1 )66 H 
( 1 0)792 2 .( 1 )66 H ■ ( 10)220 x ( 1 0) 220.(5)924 H 
(10)495i.(1)220 h .(10)792x(10)792.(1)66h 
.(l)66 H .(10)792 2 _ 

(10)792 2 .(l)66 H .(10)220x(10)220.(5)924 H 
.(l)220 H .(10)495i _ 

(10)495 1 .(1)220 h .(10)792x(10)792.(1)66h. 
.(1)66 h .(10)792x(10)792.(1)220 h .(10)495i 

Down-Type Quark Mass-Term Diagrams 



Dim 4: 


U33: 


Dim 5: 


U23: 




U32: 


Dim 6: 


U13: 




UjIi 




U22: 


Dim 7: 


U12: 




U21: 


Dim 8: 


Ull: 



(10)2203 
(10)792 2 

(10)220x (10)220.(5)924 H .(10)220 3 
(10)792x (T0)792.(l)220 H .(10)495i 
(10)220x (10)220.(1)66 H .(10)792 2 
(10)220x (10)220.(5)924 H .(10)220x (10)220 

(10)220x (T0)220.(l)66 H .(10)792x (10)792 

(10)220x (10)220.(5)924 H .(10)220x (T0)220 



Dim 5: 



Dim 6: 



Dim 7: 



Dim 8: 



D32 
D33 
D31 
D22 
D23 
D12 
D21 
D13 
Dll 



( 1 0)2203 . (5)924 H . (5 ) 220 x (5 )220.( 1 )66 H . (5)792 2 

(10)220 3 .(5)924h.(5)220x(5)220.(1)66h.(5)7923 

(10)220 3 .(5)924h.(5)220x(5)220.(1)66h.(5)792x(5)792.(1)220h.(5)66 1 

(10)792 2 .(1)66 h .(10)220x(10)220.(5)924 h .(5)220x(5)220.(1)66 h .(5)792 2 

g0)7^.(l)66 H .(10)220x(10)2m(5>24H.(5^ 

(10)495i.(1)220h.(10)792x(10)792.(1)66h.(10)220x(10)220.(5)924h.(5)220x(5)220.(1)66h.(5)792 2 

(10)792 2 .(1)66 h .(10)220x(10)220.(5)924 h .(5)220x(5)220.(1)66 h .(5)792x(5)792.(1)220 h .(5)66i 

(10)495 1 .(1)220h.(10)792x(10)792.(1)66 h .(10)220x(10)220.(5)924 h .(5)220x(5)220.(1)66 h .(5)7923 
(10)495 1 .(l)220 H .(10)792x(10)792.(l)66 H .(10)220x(10)220.(5)924 H .(5)220x (5)220 
.(l)66 H .(5)792x(5)792.(l)220 H .(5)66i 



The convention is adopted that the left-handed fields appear on the left and the left-handed conjugate fields appear on 
the right. 

(10)220 3 (10)220 3 
*> ■ ■* 

U33: i (5)924 H 



For convenience we introduce the following short-hand notation to describe each of these diagrams: 

U33: (10)220^.(5)924 H .(10)220l, D33 : (10)220^.(5)924 H . (5)220 x (5)220.(1)66 H . (5)7927, (10) 

The leading order term for U33 is seen to have dim-4, while that for D33 has dim-5, due to the (1)66h SU(5) Higgs 
singlet insertion resulting in one extra external Higgs field. 

The full sets of leading order up- and down-quark diagrams for each matrix element is presented in Table I, while 
those for the Dirac and Majorana neutrino diagrams are listed in Table II. It is rather remarkable that only one diagram 
for each matrix element appears at leading order for all four mass matrices. 



(10)220 3 (5)220 (5)220 (5)792 3 

> ? « X > ■ « 

D33: J ! (9) 



| (5)924 H A (1)66 h 



MASS MATRICES AND MIXINGS 

Given the leading-order diagrams for each matrix element in Tables 1 and 2, we can then construct the quark and 
lepton mass matrices as follows. To each diagram corresponds a coupling constant or prefactor, h"j, hfj, hf" or hy n for 
the /jth element of the appropriate mass matrix, which is assumed to be of order one at the SU(12) unification scale, 
as naturalness predicts. Every SU(5) Higgs singlet insertion in higher-order tree diagrams introduces one power of 



Table 2. Leading order Dirac and Majorana neutrino diagrams for each mass matrix element. 



Dirac-Neutrino Mass-Term Diagrams 

Dim 4: DN23: (5)7^.(5)924 H .(1)12^ 

DN33: (5)792 3 .(5)924 H .(1)12 3 _ _ 

Dim 5: DN13: (5)66i.(1)220 h .(5)792x(5)792.(5)924h.(1)12 3 

DN22: (5)792^. ( 1 )66 H ■ (5 ) 220 x (5 ) 220.(5)924 H .( 1 )220j 

DN32: (5)792 3 .(l)66 H .(5)220x (5)220 .(5)924 H .(1)220 2 

Dim 6: DN12: (5)66i.(l)220 H .(5)792x(5)792.(l)66 H .(5)220x(5)220.(5)924 H .(l)220 2 

DN21: (5)792 2 .(1)66 h .(5)220x(5)220.(5)924h.(1)220x(1)220.(1)66 h .(1)792 1 

DN31: (5)792 3 .(1)66 h .(5)220x(5)220.(5)924 h .(1)220x(1)220.(1)66h.(1)792i 

Dim 7: DN11: (5)66i.(1)220h.(5)792x(5)792.(1)66 h .(5)220x(5)220.(5)924 h .(1)220x(1)220.(1)66h.(1)792i 

Majorana-Neutrino Mass-Term Diagrams 

(l)7927.(l)66 H -(l)7927 
(1)12 3 .(1)66 H .(1)12 3 

(1)792!. (l)66 H .(l)792x(l)792.(l)66 H .(l)220 2 

(l)220 2 .(l)66 H .(l)792x(l)792.(l)66 H .(l)792i _ 

(l)792i.(l)66 H .(l)792x(l)792.(l)66 H .(l)220x(l)220.(l)66 H .(l)12 3 
(l)12 3 .(l)66 H .(l)220x(l)220.(l)66 H .(l)792x(l)792.(l)66 H .(l)792i 
(l)220 2 .(l)66 H .(l)792x(l)792.(l)66 H .(l)792x(l)792.(l)66 H .(l)220 2 
(l)220 2 .(l)66 H .(l)792x(l)792.(l)66 H .(l)792x(l)792.(l)66 H .(l)220x(l)220.(l)66 H .(l)12 3 
(1)12 3 .(1)66h.(1)220x(1)220.(1)66h.(1)792x(1)792.(1)66h.(1)792x(1)792.(1)66h.(1)220 2 



Dim 4: 


MN11: 




MN33: 


Dim 5: 


MN12: 




MN21: 


Dim 6: 


MN13: 




MN31: 




MN22: 


Dim 7: 


MN23: 




MN32: 



£ = A^su(5)/-^su(i2) ~ 1/50 through the appearance of the ratio of the singlet Higgs VEV to the mass of the conjugate 
fermion fields after the latter are integrated out. Finally as a result of the electroweak spontaneous symmetry breaking, 
the 924h acquires a weak scale VEV, v. Hence for the two quark diagrams illustrated, the matrix element contributions 
are 

U33 : «j 3 v t[t c L , D33 : h d 33 ev b T L b L L . (11) 

Note that as a result of the chiral SU(5) irrep structure, the lowest order tree diagram contribution to the 33 element 
of the charged lepton mass matrix is just the reflection of the diagram for the down quark 33 mass matrix element 
about the center of the diagram. Thus its 33 matrix element contribution is just the transpose of D33. More generally, 
the prefactors are related by h\- = h d t . By the same reasoning, it is clear that the up quark mass matrix elements are 
symmetric under interchange of i and j. 

From Table 1 we then see that the two quark and charged lepton mass matrices are given by 

h" l3 e 2 ' 

Mr = \ h'l 2 e- h'i 2 e- h u 23 e 

h 33 

h d 3 e 3 " 

Mi, = I /iV' h'Ue- h d 23 e 2 \ (12) 

h* 3 e 

h\ 3 e 2 

M L = \ h f 2l e 3 h ( 22 e 2 h[ 3 e |v = M^. 

h e 33 t 

While the up-quark matrix is symmetric, the down-quark and charged-lepton mass matrices are doubly lopsided in that 
the terms with h d 3 and h! 32 are suppressed by one extra power of £ compared with the h d 2 and h l 23 terms, respectively. 
For M D , for example, this implies that a larger right-handed rotation than left-handed rotation is needed to bring the 
down quark matrix into diagonal form, while the opposite is true for Ml. 



! n fc 


h u 12 e 3 


! 12 fc 


h\ 2 e 2 


l 13 b 


h 23 e 


,d c 3 
-21 fc 

i x e 2 


hd f\ 

h 22 e 
h d 2 e 


! 3i£ 3 


h{ 2 £ 3 
h 32 £ 2 



With the heavy right-handed neutrinos assigned to SU(5) singlets in (7), the resulting Dirac and Majorana neutrino 
33 mass matrix elements receive the following dim-4 contributions as seen from Table 2: 



DN33 : /$vv 3 lM 



*3L' 



MN33 : Kg\ R Nl L Nl L . 



(13) 



Here A R represents the right-handed mass scale, typically of <^(10 14 ) GeV, whereas the SU(5) SUSY GUT scale is 
2 x 10 16 GeV to obtain gauge coupling unification. Again, a factor of e enters for every singlet Higgs insertion in 
higher order diagrams. The two neutrino mass matrices can then be read off from Table 2, and we find 



M DN = 



M MN = 



'h dn x e 3 

"21 fc 

h™e 2 







h 22 £ 




hi»e 


hil) 


h n £ 


h^s 2 


h™e 2 


h^e 3 


h™e 3 


h 33 



(14) 



A*. 



where Mqn is also double lopsided, while Mmn is complex symmetric as usual. The symmetric light-neutrino mass 
matrix is obtained via the Type I seesaw mechanism: 



M v = -M DN M M ^Ml N . 
Keeping only the leading-order terms in e for each matrix element, we find 



(15) 
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(16) 



which does not involve the pref actors /if", ^fi > 



'21' "31- 



/z™ and A™. 



The light-neutrino mass matrix exhibits a much milder hierarchy compared to the up-type and down-type mass 
matrices, as can be seen from the pattern of powers of £. A mild or flat hierarchy of M v is conducive to obtaining 
large mixing angles and similar light neutrino masses. Furthermore, one observes that the light neutrino mass matrix 
obtained via the seesaw mechanism involves the doubly lopsided Dirac neutrino mass matrix twice. The lopsided 
feature of Mdn is such as to require a large left-handed rotation to bring M v into diagonal form. 



NUMERICAL RESULTS 

From the above up and down quark, charged lepton and light neutrino mass matrices, one can diagonalize the cor- 
responding Hermitian matrices, MM\ in the usual manner to obtain the mass eigenvalues and the unitary transfor- 
mations, U, effecting the diagonalizations. The Cabbibo-Kobayashi-Maskawa (CKM) quark mixing matrix and the 
corresponding Pontecorvo-Maki-Nakagawa-Sakata (PMNS) lepton mixing matrix then follow as 

V C KM = UlU D Vpmns = uIu v (17) 

To obtain numerical results for the model predictions, we evaluate the mass matrices at the top quark mass scale and 
use just real prefactors to avoid too many fit parameters for good fit convergence. There are 6 prefactors each for the 
symmetric up quark and Majorana matrices, 9 each for the lopsided down quark and Dirac neutrino matrices, but 5 of 
them do not appear in the light neutrino mass matrix, making a total of 25 parameters. In addition, we have one for the 
right-handed neutrino scale, Ar plus a value for e which we fix at e = 1/6. 5 2 = 0.0237 again for good fit convergence, 
for a grand total of 26 adjustable fit parameters. To avoid correlated data parameters, we make use of the 9 quark and 
charged lepton masses plus the 3 neutrino Am 2 's and the 18 CKM and PMNS mixing parameters taken to be real, for 
a total of 30 data points. We refer the reader to our published paper [13] for full details of the fitting procedure, where 



we have included the latest best value for the reactor neutrino mixing angle, 613. There can be found a table giving the 
phenomenological mass and mixing data entering the fit, as well as the theoretical mass and mixing results obtained 
from the fitting procedure. 

The best fit was obtained with a normal neutrino mass hierarchy with Ar = 7.4 x 10 14 GeV and the following mass 
matrices: 

/-Lie 4 7.1e 3 5.6e 2 \ /-6.3e 4 8.0e 3 -1.9e 3 \ 

M v = ( 7.1e 3 -6.2e 2 -O.lOe )v, M D = ( -4.5e 3 0.38e 2 -1.3e 2 )v =M[, 



(18) 
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5.6e 2 
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-O.lOe 


-0.95 


h^e 3 


0.21e 2 


-2.7e\ 


h d 2 1e 2 


-0.28e 


-0.15 h 




2.1e 


0.086 / 


-81e 2 


-4.3e 


2.4e \ 


-4.3e 


-0.25 


0.28 v ; 


2.4e 


0.28 


-1.1/ 



M DN = \ h d "e l -0.28e -0.15 v, M MN = 



M v = -4.3e -0.25 0.28 w 2 /A R . 




Note that all prefactors except three in the above matrices are within a factor of 0.1 to 10 of unity for this best fit. The 
five independent prefactors, hf", hf", hf{, h™" and h™", do not influence the fit and remain undetermined as noted 
earlier. For this best fit we find the neutrino mass values 

mi=0, m 2 = 8.65, m 3 =49.7meV; M x = 1.67 x 10 12 , M 2 = 6.85 x 10 13 , M3 = 5.30 x 10 14 GeV. (19) 

In addition, the best fit favors 8cp — % for the leptonic CP Dirac phase. The value of e used then implies that the 
SU(12) GUT scale is about M SU ( 5 j/e = 8.4 x 10 17 GeV, just below the reduced Planck scale, where we have used 
2 x 10 16 GeV for the SU(5) unification scale. 

All remaining mass and mixing parameters are fit quite well by the model; however, since Ml is just the transpose of 
Mo in leading order in e, the Georgi-Jarlskog relations [17] are not satisfied for the down quarks and charged leptons. 
We have checked that the addition of an adjoint 143h Higgs field whose VEV points in the B — L direction contributes 
to M D and M L at one higher order of e, so that the M L = M^ relation is broken, and more accurate values can be 
obtained for the down quark and charged lepton mass eigenvalues. 



SUMMARY 

A unified SU(12) SUSY GUT model was obtained by a brute force computer scan over many SU(12) anomaly-free 
sets of irreps containing 3 SU(5) chiral families under the assumption that the symmetry is broken in stages from 
SU(12) — > SU(5) — > SM. In doing so, looping over all SU(12) fermion and Higgs assignments was performed with 
good fits to the input data required. For this purpose an effective theory approach was used to determine the leading 
order tree-level diagrams for the dim-(4 + n) matrix elements in powers of e", where £ is the ratio of the SU(5) to the 
SU(12) scale. The best fit was obtained by requiring all prefactors to be but the large number of them implies 

just a few predictions. With no discrete flavor symmetry adopted, problems with breaking by gravity, domain walls 
and explanation of its origin can be avoided [18, 19, 20]. On the contrary, with such a large SU(iV) gauge group, a host 
of heavy fermions is predicted which are integrated out at the SU(5) scale. 

The SU(12) model considered is just one of many possibilities (including other assignments and larger SU(7V) 
groups), but its features were among the most attractive found: Each SU(5) family supermultiplet can be assigned to a 
different SU(12) multiplet in the anomaly-free set. In the model considered, only one diagram appears for each matrix 
element for all 5 mass matrices, but some additional contribution is needed to obtain the Georgi-Jarlskog relations. 

Among some distracting features we point out the prefactors are determined at the top quark scale. They should be 
run to the SU(5) unification scale to test their naturalness. The fit considers only real prefactors, so CP violation is 
not accommodated, but the fit preferred Sep = % over 8cp — for the leptonic CP phase. The complete breaking of 
SU(12) — > SU(5) while preserving supersymmetry needs to be worked out in more detail and is under further study. 
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